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^) ' We calculate the two-loop diagrams that give a non-zero contribution to the 

js^ ■ electric dipole moment di of a charged lepton / due to possible Majorana masses of 

tH- , neutrinos. Using the example with one generation of the Standard Model leptons 

^D I and two heavy right-handed neutrinos, we demonstrate that the non- vanishing result 

for di first appears in order 0{mi7n'^G'p), where nii, is the mass of the light neutrino 

Mh, and the see-saw type relation is imposed. This effect is beyond the reach of presently 

Q^' planned experiments. 
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1 Introduction 

Recent discovery of CP violation in the neutral 5-meson system [1, 2] is in perfect accord 
with the CP violation observed in K mesons [3]. Both results are explained within the 
minimal model of CP violation, known as the Kobayashi-Maskawa model [4]. It links 
CP-violation in K and B mesons to a single CP-violating invariant of the Kobayashi- 
Maskawa matrix in the quark sector, Im(\4;,Vj^VcrfV^|,) ~ 3 x 10~^. This combination, as 
well as ^QCD, are the only sources of CP violation in the Standard Model. 

An independent piece of experimental information about CP violation comes from the 
searches of the electric dipole moments (EDMs) of neutrons [5] and heavy atoms [6, 7]. 
It is well known that the Kobayashi-Maskawa model predicts extremely small EDMs. 
Indeed, the necessity of four electroweak vertices requires from any diagram capable of 
inducing an EDM of a quark to have at least two loops. Moreover, it turns out that 
all EDMs or color EDMs of quarks vanish exactly at the two-loop level [8], and only 
three-loop diagrams survive [9, 10], producing a tiny number of order lO^^^e cm. 

The only relevant operator that is not zero at the two-loop order is the Weinberg operator 
[11], but its numerical value turns out to be also extremely small. Possible enhancement 
comes from the large distance effects, that could lead to a KM-generated EDM of the 
neutron of order 10~^^e cm [12], which is still six to seven orders of magnitude smaller 
than the current experimental limit. The KM phase in the quark sector induces the EDM 
of a lepton via a diagram with a closed quark loop, but a non-vanishing result appears 
first at a four-loop level [13] and therefore is even more suppressed. 

The suppressed values of EDMs from the KM model together with enormous accuracy 
of EDM experiments produce stringent constraints on possible flavour-diagonal sources 
of CP violation, notably on CP-odd combinations of the soft-breaking parameters in the 
super symmetric extensions of the Standard Model (SM). 

On the other hand, the KM phase cannot be the only source of CP violation in nature. 
Dynamical generation of the baryon asymmetry of the Universe (BAU) requires presence 
of an additional source(s) of CP violation in nature. One of the most appealing scenarios 
for BAU is leptogenesis [14] where a non-zero lepton number is generated from the out-of- 
equilibrium decay of heavy right-handed Majorana neutrinos with CP violation coming 
from the Yukawa sector of the theory. Subsequent sphaleron transitions [15] transform 
half of the initial lepton asymmetry into BAU. 

The attractiveness of leptogenesis is in its simplicity, relative freedom from the low-energy 
constraints, and in recent experimental results that indicate large mixing angles in the 
lepton sector [16]. 

The complex phases in the lepton mixing suggested by leptogenesis at some level will 
induce de-, the electron EDM, as well as EDMs of other charged leptons. Various calcula- 
tions of df. and rf^ were performed in the supersymmetric case under certain assumption 
about the soft-breaking parameters [19], while a non-supersymmetric case remains poorly 
explored. 

We note that manifestations of Majorana phases in CP-violating phenomena has been 
explored in [20]. The EDM is an additional observable arising due to those phases. Com- 
pared with the KM model, where three generations and at least four loops are necessary 



to generate a lepton EDM, the addition of Majorana phases enhances the effect. However, 
as we will see, it remains very small. 

In this paper, we present a systematic analysis of the two-loop diagrams that lead to the 
ED Ms of charged leptons through interactions with two additional heavy neutrinos, in a 
one-generation model suggested by a minimal leptogenesis scenario. These diagrams can 
be divided into two major classes. The first class where the lepton number is conserved 
along the fermion line vanishes identically, as the same arguments that lead to dq = at 
two loops [8] apply (see Fig. 1). 





(a) (b) 

Figure 1: Examples of two- loop diagrams that give zero electric dipole moment of the 
electron. The direction of the arrows remains the same on all fermion propagators and 
the external photon can be attached to any charged line. 



The second class of diagrams (Fig. 2) has no analogues in the quark sector: due to 
the existence of the Majorana component of the neutrino mass matrix, a combination of 
AL = 1 and AL = — 1 transitions is possible. The existence of this additional class of 
diagrams for leptons was pointed out in Ref. [17]. It leads to a non-zero electron EDM 
{de) even at the two-loop level, but up to now there has not been any serious attempt to 
calculate the size of de induced by this class of diagrams. (An estimate of d^ due to CP 
violation in the sector of heavy Majorana neutrino was presented in Ref. [18]. However, 
this result is not satisfactory, as it suggests that de does not vanish in the limit of infinitely 
heavy right-handed neutrinos.) 

The purpose of this work is to give a detailed calculation of de due to the CP violation 
in the lepton sector in the presence of Majorana masses for neutrinos and obtain a pre- 
diction for the electron EDM in the see-saw model of the neutrino mass sector. A similar 
mechanism can of course also lead to the EDM of the muon. 

We find that two-loop diagrams give non-vanishing results for the de- However, EDMs 
are very much suppressed by the smallness of the neutrino mass as they first appear 
in OimetnlGy) order, if the see-saw type relation between rriy, Dirac mass and heavy 
Majorana mass is imposed. Thus, for any phenomenologically motivated choice of rrii, the 
result for de (and d^) turns out to be much smaller than current or projected experimental 
sensitivity to EDMs. We also note that the fine-tuning in the neutrino mass sector may 
lead to up to ten orders of magnitude enhancement in de-, which is still not enough to 
bring it within the experimental reach in the near future. 
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Figure 2: Contributions to the electron EDM in a model with Majorana masses of neu- 
trinos. /i,2 denote all possible neutrinos (see text). Crosses denote insertions of lepton- 
number violating mass parameters. Note that the direction of the internal electron line 
is opposite to the external ones. 

2 Description of the model 



We take one standard model generation: , cr, and two singlet heavy neutrinos 

Ni^2- The latter do not participate in electroweak interactions; in particular, the charged 
current sector is described by the Lagrangian 



£,, = -^ {uLVteL + (H.c. 
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The mass sector Lagrangian for fermions is 

— ^M 

+ 

+mi [e"^'NiUL + c-'^^vlNi 
+m2 (e^^^iVsz/L + c-^^'vlN^) 



2 ^NINi + NiNI) + '-^ [n^N2 + N2N^ 
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Here ip'^ = 'yoCip*; Mi^2 and mi^2 are defined in terms of real positive Yukawa couplings 
yi^2 and the electroweak vacuum expectation value v, 
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We use the freedom of phase choice for ul and cr^l to redefine 



z/i^e-'^V^. 



(3) 



(4) 



We see that there is only one physical CP violating phase 77 = 0i — 02- 

Before we explore the physical manifestation of rj, we determine the mass eigenstates of 
neutrinos. We will use them to compute the EDM of the electron. We use the identity 






-[v^N + N^v': 
^ 'NvL + ylN" 
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to rewrite the neutrino mass matrix as 
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We use Ml 2 ^ wii 2 to approximately diagonalize the neutrino mass matrix. We define 
neutrino mass eigenstates: hght z/ and heavy ni 2, 



Ni \=V\ ni 

N2 J \n2 

where V is the matrix diagonahzing Ai. Approximately, 
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where we denote small corrections by o„ = (9 (m ( ^ ) ) . Two things are worth noting: 
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Heavy neutrinos ni^2 participate in charged current interactions described by eq. (1), 
through their presence in z/^, eq. (7): 
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The CP violating phase enters through the complex mass of v and through the 
charged current interaction of n\. 



The Majorana mass of the light (active) neutrino is given by the following relation, 
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Since experiments with light neutrinos give information about rriiy, it is convenient to keep 
this parameter fixed, while allowing 7711^2, ^1,2 and rj to vary. In doing this, we would 
have to distinguish two possibilities. 



See-saw relation, when both contributions in (10) are on the order or smaller than 



• Cancellation of two terms in m,y which we would term as a "fine-tuned case", al- 
though it could be a result of some symmetry that suppresses the determinant of the 
mass matrix (6). For this cancellation to happen, one has to have m\/Mi ~ m\/M2 
and phase 7] close to 7r/2. 

We would like to note that theoretically one could have a fine-tuned case even for mi,2 ~ 
Ml 2- This, however, would also result in large admixtures of heavy rii and n2 neutrinos 
in the original left-handed neutrino. In a more realistic model with three active neutrino 
species, such admixtures could lead to a non-universality of charged currents, modification 
of the invisible Z decay width, etc. With the accuracy of electroweak data, we impose 
the bound on the size of the mixing angle: 

— , — <0f-V (11) 

Ml' M2 ~ uoy ^ ' 

3 Computation of the electron EDM at two loops 

To define the electron EDM, consider the general matrix element of an electromagnetic 
current between spin 1/2 fermions. 



+75 

with 1 = 0^. The EDM de of an electron of charge e is given by 



2m m 



u,ip + q), (12) 



de = -TZ^G^iO). (13) 

Zrrie 

The most notable property of all diagrams in Figs. 1 and 2, potentially contributing to the 
CP-odd amplitude, is the effective antisymmetrization over the neutrino propagation, /i 
and /2. This is because the CP-odd part of a diagram with the selection of eigenstates /i 
and /2 is always opposite to the CP-odd part of the diagram with interchanged flavours. As 
a result of this antisymmetrization, in the expansion over small q the flrst class of diagrams. 
Fig. la and lb, turns out to be proportional to the cube of the photon momentum O(g^), 
while 6*2(0) vanishes identically. A more detailed explanation of the cancellation of EDMs 
can be found in Ref. [8, 13, 11]. 

Within the model we are considering, the non-vanishing diagrams contributing to de 
are shown in Fig. 2. Their evaluation is simplifled if we assume the most natural mass 
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hierarchy, ^TI "^ ^^ < Mw < Mi,2- We use the notation M = M±M and AM = 
M2 — Ml and assume |AM| ^ M for calculational convenience. 

First, we treat the case when /i 2 are both heavy (some details of derivations of these 
results are given in the Appendix), 

. w, , X fGF\^ AMm^m^ /16 364 112 .\ , ,, 

A4 heavy-heavy = e ^^ ^e^^^TT^ + ^-vrM sin 2r7. 14 

^^ ^ ^' V167r2y ^ M M M V3e 9 27 / ' ^ > 



Next, we obtain the contribution when one of /» is the hght state z/, 

A4(hea,vy-light) = e (^) rn,—^^M-<-M^ (-- + -) sm2„. (15) 
The sum of both contributions is finite, 



/GfX^ AMm2m2/32, M 260 112 ,n 
e^^ me— — — i— ^ — In— + ^-ttM sin 2r/. (16) 



df, = A(ie (heavy-heavy) + A(ie(heavy-hght) 

f Gf \^ AMm^m^ /32 , M 

r mp — In 

V167r2y ' M M M V 3 Mw 9 27 

Our answer shows that the electron EDM first appears at 0{rneTn\'m\M~'^G']p) level, which 
is essentially the same as 0{'me'mlG'^) if the see-saw relations 771^2/^ ~ 0{my) are in 
place. 

In a slightly different model, with two generations of active neutrinos with masses niyi 
and 711^2 there is an additional possibility of /i and /2 both being light. It is easy to see, 
however, that in the case of (light-light) diagram the contribution to d^ is parametrically 
much smaller than Eq. (16), being suppressed by four powers of light neutrino mass 
myimj^2{f^'ii ~ ^u2)- Thus, the generalization to the case of three generations would not 
bring any parametric change to the size of d^- Fixing rriu, we can also allow mi 2 and Mi 2 
to vary. If Mi^2 becomes smaller than the electroweak scale, de acquires additional sup- 
pression by M'^ /My^. This proves that (16) with the chosen hierarchy of scales represents 
the largest possible contribution to de at fixed rriy. 



4 Discussion and Conclusions 

What is the largest numerical value of the answer (16)? For the see-saw type relation, 
in the absence of substantial cancellations between the two terms in (10), we can use the 
following inequality 



2 2 
mf mn 

sm2r7 

M M ' 



<mt. (17) 



This leads to an extremely tight theoretical bound on the possible size of the EDM, 



fGF\'.\AM\f^^^^^^^ M 



Mel ^ e mem,, 10.7 x In h 12.1 , see-saw case. (18) 

Using |AM| ~ M ~ 10^^ GeV, and the Particle Data Group book [21] bound on the mass 
of the electron neutrino, iriij < 3 eV, we arrive at our final numerical result for the see-saw 



Mel < 1.5 X 10 ^^e cm, see-saw case. (19) 

Not only is this number much smaller than the most optimistic accuracy of future electron 
EDM searches [22], but also it is significantly smaller than EDMs induced by the KM phase 
from the quark sector. 



The fine-tuned case deserves special consideration, as the numerical answer for de can be 
significantly larger than (19). Introducing another angle 5 = tt/2 — rj, which has to be 
small for the fine-tuning of rrij, to happen, we rewrite the light neutrino mass squared as 

^' = ^2 (K - mlf + 4mlml6') . (20) 
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On the other hand, the answer for EDM contains a factor — pr^ sin(2r7) = 2^^6, which 
in view of Eq. (20) can be limited as 






< ^^^^- ^ (0-1 X 500 GeV)m,. (21) 



In this expression, we have used (11) and the upper bound on mi 2 taken to be ~ 500 
GeV. The latter follows from the condition that Yukawa sector remains perturbative, i.e. 
1/1,2 ^ 0(1 — 10). Comparing (17) and (21), we observe that the fine-tuned case may lead 
to up to 10^° enhancement relative to the see-saw case and the electron EDM may reach 
a value of 

de < 10^^^ e cm, fine-tuned case. (22) 

This value is still far from the existing or projected experimental accuracy but is much 
larger than the electron EDM induced by the KM phase. 

In some models the effects of the CP-odd electron-nucleon interaction, CsNNei'-f^e may 
dominate over the electron EDM contribution [23] in the atomic (molecular) EDM. An 
analysis of possible two-loop diagrams for Cs shows that in the model with CP violation 
coming from Majorana neutrinos this is not the case, and the contribution of Cs{r)) to 
atomic EDM is smaller than (19) for the see-saw case and (22) for the fine-tuned case. 

To summarize, we have calculated the contributions of the two-loop diagrams to the EDM 
of the electron in non-sup ersymmetric models with CP violation in the lepton sector and 
with Majorana masses for neutrinos. We notice that the non-zero result for de can be 
obtained in a rather minimalistic way: with one generation of SM leptons and two right- 
handed neutrinos. In terms of the right-handed mass M and Dirac mass rriD, the non-zero 
result appears in order 0{me'mj^M~'^G'jp), which is the largest possible result given the 
symmetries of the model. If the smallness of the light neutrino mass is achieved via 
usual see-saw relation without any fine-tuning, this answer is equivalent to 0(mem^G|.) 
and numerically is extremely small. Numerical smallness of this result stems from the 
smallness of rriiy, and the parametric suppression of EDM is very similar to the suppression 
of e.g. the SM amplitude for fi —^ e'j decay. In the fine-tuned case, when the smallness 
of the light neutrino mass is achieved via near perfect cancellation of two contributions, 
the result for de may become larger by many orders of magnitude, but still smaller than 
about 10~^^e cm, and therefore much smaller than the sensitivity of any EDM experiment 
in the foreseeable future. 

We conclude that the KM-type models with Majorana neutrino masses do not have any 
impact on the EDM searches. Therefore, the only options of searching for CP-violation 
in the neutrino sector suggested by the leptogenesis scenario are the CP asymmetries in 
neutrino oscillations and not EDMs. Conversely, possible positive results from the future 
electron EDM searches could be an indication of leptogenesis combined with the soft 
supersymmetry breaking. 
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A Details of the calculation of two-loop diagrams 

Here we give intermediate results for the sum of the three diagrams in Fig. 2, before we 
add and simplify contributions from ni^2- 

For the heavy-heavy case we get 

ie^ 1 / 17111712 \ ^ / 2?n \ 

A4(heavy-heavy) = -^^^^(I^;^ (^^^^^j (e ''- complex conjugate) 

..^ Mmj ,, A, fl 91 7 2' 



M^ V6e 72 54 



f Gf \^ AMm^m^ /I6 364 112 ., , • . 
' > mec ^^M-^M + TT^ sin2?7. 



Vl67r2y M3 V3e 9 27 

(23) 



For the heavy-light case we find 
A(ie (heavy-light) = 



te 1 mz 



^ e^^V f-^' ^il + ^ / f-^e-2^^ M2] - complex conjug; 
Mf ■' \ M2 J Mi \ Ml J 



(24) 



where the function / arises from the evaluation of the diagrams without coupling con- 
stants, 

....4ef ^ . Mf 5 13, Mf l,2^f 31 7r2\ 

/(m^. Mi) = m^MiMw' n -^ + — n —^ + - n^ -^ + — • 

^^ ' ^ ^ \^ 12e M^ 16e 72 M^ 8 M^ A8 36 J 

(25) 

Here rrii, and Mi denote a light and a heavy neutrino mass insertion. Working to first 
order in the heavy neutrino mass splitting AM, we get 

ie^ 1 ml 



A(ie(heavy-light) 



2m,,s4 (167r2)2 Ml 



m?m2 AM,^4 / 16 104 _ , M^ \ ^ 

— ^—^i M^^' + + 16 n^ sn2r7 

M3 16 '^ V 3e 9 Mlrj ' 

f Gf \^ AMm?mL . / 16 104 , M^ \ 

(26) 

We see that in the sum of contributions in (23) and (26), the divergences cancel and we 
obtain the final result, eq. (16). 
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B Asymptotic Expansions 

In this appendix, we show how the method of asymptotic expansions [24] is used to 
calculate the contributing diagrams for the electric dipole moment. 

Two-loop Feynman diagrams containing more than one mass scale are difficult to solve 
analytically. Analytic expansions in small parameters, such as the ratio of lepton and 
weak boson masses, are used to reduce the calculations to integrals involving only one 
mass scale. The results are integrals which have analytic solutions. As an example of this 
technique, we use Feynman diagrams from Fig. 2 containing one light neutrino, /i, of 
mass rriiy and one heavy neutrino, /2, of mass M. 

The momentum assignments for the calculation are illustrated in Fig. 3. Conservation of 
momentum at each vertex requires pi = p + ki, p2 = p + ki + k2 and ps = p + k2- 




P Pi P2 P3 P 

Figure 3: Momentum assignment for the EDM calculation. 

In order to perform the expansions, all mass scales must be identified. For our calculation, 
the mass scales involved are M ^ My/ ^ nif. ^ niy. 

Once the mass scales are determined, the integral volumes are divided into regions so that 
the momentum fiow through the internal lines are on the order of one of the mass scales. 
With one light and one heavy neutrino, there are four scenarios to be considered: 

1. ki,k2<^M,Mw 

2. k2 ~ M, fci < M, Mw 

3. ki,k2 ~ M 

4. ki, k2 ~ M, /ci + A;2 ~ nie 

Within each region, the appropriate propagators can be expanded in a Taylor series to 
reduce the number of mass scales in the integrals. 



B.l ki,k2<^M,Mw 

In this momentum region, the heavy neutrino and the W boson propagators can be 
expanded in a Taylor series. For example, the heavy neutrino propagator can be written 

-^-^y(J^Y (27) 
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since k <^ M. 

The Taylor expansion reduces the diagram in Fig. 3 to that of Fig. 4. The two-loop 
diagram is reduced to the product of two one-loop diagrams. Each one-loop diagram has 
one mass scale and has an analytic solution. 




Figure 4: Reduction of a two-loop diagram to two one-loop diagrams for ki, k2 <^ M, M]j\r. 
The solid lines represent massive propagators. 



B.2 A;2 - M, A;i < M, MvK 

When one momentum is large and the other is small, the propagators are expanded in all 
momenta except for the large momentum. This reduces the calculation to the diagrams 
in Fig. 5, a product of two one-loop diagrams. Clearly, the first diagram has only one 
mass scale involved. In the second diagram, partial fractions are used to further separate 
the mass scales. In particular, we use the identity 



1 1 



fc2 + M2 A;2 + M^ M^ - M2 \k^ + M2 P + M^ 
The resulting diagrams can be computed analytically. 



M 



:28i 



V_y ^Taaat^ 



niv 



Figure 5: Reduction of a two-loop diagram to two one-loop diagrams for k2 ~ M, ki <^ 
M,Mw- The dashed line represents a massless propagator. 



B.3 ki,k2-^M 

If both momenta are large, a reduction is made to the form of Fig. 6. Notice however, 
that there are still two mass scales in the diagram. Thus, we use the assumption that 
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M ^ Mw for further expansion to reduce the diagram down to a one mass scale diagram. 



Figure 6: Taylor Expansion for ki^k2 ~ M. 



B.4 /ci, A;2 ~ M, /ci + A;2 ~ me 

The last case reduces to Fig. 7 but does not contribute to the d^ in the leading order. 




Figure 7: Taylor expansion for the case ki + k2 <^ M. 

After integration over the initial ranges is performed, the contributions from all momen- 
tum regions are summed to produce eq. (15). 

Using a similar analysis, eq. (14) was computed for the diagrams with two heavy neutri- 
nos. 
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